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Scalar QED is studied with higher order derivatives for the scalar field kinetic energy. A local 
potential is generated for the gauge field due to the covariant derivatives and the vacuum with 
non-vanishing expectation value for the scalar field and the vector potential is constructed in the 
leading order saddle point expansion. This vacuum breaks the global gauge and Lorentz symmetry 
spontaneously. The unitarity of time evolution is assured in the physical, positive norm subspace and 
the linearized equations of motion are calculated. Goldstone theorem always keeps the radiation field 
massless. A particular model is constructed where the the full set of standard Maxwell equations is 
recovered on the tree level thereby relegating the effects of broken Lorentz symmetry to the level of 
Q_i radiative corrections. 
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I. INTRODUCTION 



We know no exact equation of motion in physics, all laws are inferred by ignoring some loosely attached part of the 
system considered. As a result, the equations of motion should be tested for the stability of solutions against adding 
Oh| small correction terms to the equations. Such an analysis is usually performed in the framework of the renormalization 
group [1] and perturbation expansion can be used to establish that perturbing terms with higher mass dimension (we 
use units c = H = 1) are less important at short distances. 

Nevertheless, it is an important difference whether the higher mass dimension arises from field amplitude or space- 
time derivative because the latter may modify the tree-level, normal mode structure and generate new degrees of 
freedom. In fact, a field theory with a real, single component scalar field characterized by a Lagrangian containing 
rid space-time derivatives of the field contains rid degrees of freedom. Once the new propagating degrees of freedom 
are present their interactions might well be non-negligible due to IR or UV divergences even if the coupling constants 
■ in the bare Lagrangian are weak. There is another more fundamental change generated by these terms, spontaneous 
symmetry breaking of space-time symmetries due to an inhomogeneous condensate, the subject of this work. The point 
is that the renormalization group equations usually include quantum fluctuations only. The higher order derivatives 
terms may generate new relevant operators in the IR on the tree-level which lead to a vacuum with inhomogeneous 
condensate. We do not embark on a general renormalization group study here, rather present a simple-minded analysis 
of the symmetry and the quasi-particle content of an Abelian gauge model in the leading order saddle point expansion. 

If the condensate consists of bosons with non-vanishing momentum, filling up the whole quantization volume, then 
the "wavy vacuum" breaks the space-time symmetries, in a manner similar to solids where the infinite inertia of 
the solid prevents the zero modes to restore the broken external symmetries. The result, expected from solid state 
physics, is the appearance of several branches of the dispersion relations, different elementary excitations in the theory. 
Note that if translation invariance is broken at sufficiently short length scale to remain undetectable for the class of 
observables one uses then the vacuum appears homogeneous. We shall see that in models with gauge symmetry where 
the covariant derivative is supposed to acquire non- vanishing value in the condensate the inhomogeneity of the vacuum 
may be gauged away and we find a homogeneous condensate which simplifies the model enormously. The result is 
some kind of extension of Higgs-mechanism where the non-vanishing expectation value for the gauge field breaks 
Lorentz symmetry. The resulting Goldstone modes remain in the gauge field sector and protect some components of 
the gauge field against mass generation. 

The model studied in this work is scalar QED where higher order (covariant) derivative terms are introduced for the 
charged scalar field. The higher order terms of this model can be imagined either as smooth cutoff in defining an UV 
finite theory or as originating from the elimination of some heavy particle and approximating the self energy of a scalar 
charged particle by a polynomial of finite order in the momentum. Goldstone theorem protects the electromagnetic 
field against becoming massive, Maxwell equations are recovered in the linearized equation of motion, rendering the 
Lorentz symmetry breaking effects to radiative corrections. The rather technical problem of proving unitarity of the 
model within the physical, positive norm subspace is solved within perturbation expansion by assuring real energy 
spectrum for normal modes and preserving the physical subspace, consisting of states of positive norm during the 



2 



time evolution. 

The dynamical breakdown of space-time symmetries by higher order derivatives has already been studied in two [2] , 
three [3] and four [4, 5] dimensional Euclidean models where periodically modulated condensate has been observed 
and several particle modes have been found corresponding to a single quantum field [5]. The present work can be 
considered as continuation of such inquiries for models defined in Minkowski space-time and equipped with gauge 
symmetry. The spontaneous breakdown of relativistic symmetries has been considered within the scheme of emerging 
photons [6] and the bumblebee models [7] where an external Mexican hat potential is assumed for the vector bosons. 
Our plan is less ambitious and starts with photons as elementary particles. 

Our results can be best summarized by comparing them with the conventional Higgs-mechanism where Goldstone 
mode arising from the spontaneous breakdown of global gauge invariance appears in the gauge field which becomes 
massive. In our case the relativistic space-time symmetry is broken spontaneously as well, leaving behind three more 
Goldstone modes. Two of them are non-vanishing helicity components of the gauge field and restore the conventional, 
massless radiation field of electrodynamics. The third soft mode resides in a certain combination of the vanishing 
helicity component of the photon and the scalar field and is responsible for the preservation of the usual, long range 
Coulomb propagator for the temporal component of the gauge field. Therefore, despite the spontaneous breakdown of 
internal and external symmetries the free propagator and the normal modes of the electromagnetic field are equivalent 
with those of conventional electrodynamics. The symmetry breaking influences radiative corrections and the dynamics 
of the charged scalar field only. 

We start in this paper in section II by listing few salient features of scalar models with higher order derivatives. 
The issue of unitarity and the way it can be recovered by proving reflecting positivity in Euclidean space-time is 
discussed in Section III. Our model, scalar electrodynamics with higher order derivative for the charged scalar field 
is introduced in Section IV. The dynamics is discussed in static temporal gauge where the exceptional features of 
the time component of the gauge field can be dealt with in the easiest manner. Section V contains the construction 
of the vacuum in the leading, tree-level order of the saddle point expansion. The stability of the vacuum and the 
unitarity within the physical subspace is shown in Section VC. The particle content of the theory is defined by the 
quadratic part of the Lagrangian which is explored in Section VI. Finally, Section VII is reserved to our summary. 
An Appendix contains the details of calculating the quadratic part of the action. 

II. UNITARITY AND HIGHER ORDER DERIVATIVES 

Effective theories may or may not be unitary. In fact, the unitarity is lost when a particle, retained in an effective 
theory can lower its energy by the emission of other particles which have been eliminated in deriving the effective theory. 
Nevertheless, non-unitary effective theory remains a powerful approximation scheme when these decay processes are 
kinematically suppressed and make the life-time sufficiently long. But one would still prefer to recover the simplicity 
following from unitarity in effective theories which tend to be rather complicated. For instance processes whose energy 
remains below the mass M of the particle eliminated should reflect unitary dynamics when considered for sufficiently 
long time. Nevertheless the UV divergences and quantum anomalies of the underlying theory mix the high energy 
effects into low energy sector. The most natural way of recovering a unitary effective theory is to place the UV cutoff 
below the eliminated particle mass, A < M. But this solution is not as simple as it seems. On one hand, smooth 
cutoff allows decay processes with small but non-vanishing probability, and on the other hand, sharp cutoff leads to 
artificial non-local, acausal dynamics at the length scale A -1 which is observable in this case. 

The hallmark of effective theories is the appearance of higher order derivatives in the Lagrangian, reflecting 
momentum-dependent self energies of quasi-particles or form factors. The latter appear in vertices and have mainly 
perturbative effects. But the former are in the quadratic part of the action in the fields and modify the structure of 
quasi-particles. They are sometimes used as a kind of Pauli-Villars regulator which renders the effective theory UV 
finite [8-10]. Even though the scale of this smooth cutoff is M, the non- unitary processes are not fully suppressed. 
Once the effective theory is rendered UV finite we may consider it as an extension of the class of potentially inter- 
esting consistent, microscopic models because its UV dynamics is well defined. Motivated by the search of possible 
fundamental theories one naturally expects the complete suppression of non-physical, non-unitary processes. 

We consider in this section a model for a neutral scalar particle described by the field <j>(x). The interaction vertices 
will be kept in a momentum independent fashion only and the Lagrangian 

C=^(x)L(-D)^x)-V(cp 2 (x)) (1) 

where the real function L(p 2 ) represents the sum of the kinetic term and a momentum dependent self energy and is 
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supposed to be a polynomial of order (p 2 ) nd which assumes the form 



L{p 2 )-V(p) = Z-'\{{p 2 -m 2 n ) 1 (2) 



n=l 



where Z is real, the potential has a minimum at </> = cfi and the poles might appear in complex pairs. The role of the 
poles p 2 — to 2 can be seen more clearly by means of partial fraction decomposition [11], 

E ( 3 ) 



L{p 2 )-V'{4> 2 ) ^p 2 -m 

where z n — Z / dL{m 2 ) / dp 2 . We assumed single roots in this equation. In case a root p 2 — m 2 is of £-th order then 
the right hand side may contain terms z^/(p 2 — m. 2 ) fc with 1 < k < £. Complex roots produce complex contributions 
to the loop-integrals and lead to exponentially damped or increasing amplitudes in time and unitarity can be saved 
by a graph-by-graph modification of the theory only [12]. Another problem is seen for real roots when the kinetic 
term L(p 2 ) is a real function and it displays slope with alternating sign at its roots. Thus approximately half of the 
contributions to the kinetic energy has the wrong sign, indicating that the Hamiltonian is unbounded from below. 
This instability can be cured by introducing negative norm states [10] but the unitarity within the physical, positive 
norm subspace could not be established in a nonperturbative manner [13]. Therefore a careless truncation of the self 
energy may spoil unitarity and stability of the effective theory. 



III. UNITARITY AND REFLECTION POSITIVITY 



A proposition to preserve the desired properties was put forward by starting with an effective theory in Euclidean 
space-time [14], where it is usually derived perturbatively. The effective theory (1) should have a well defined Euclidean 
path integral representation, a condition assured by imposing the constraint L{— p^ 2 — p 2 ) > 0. The safest is to use 
lattice regularization in Euclidean space-time where higher order derivatives can be represented as higher order finite 
differences. It is easy to see in lattice regularization that we need new variables to regain the usual description for 
theories with higher order derivatives [15]. The Kolmogorov-Chapman equation expresses the group structure of the 
time evolution in the Fock-space and can be written as 



tl [0 (3 \* (1) 



J £)[^ 2 )]e- St 3-*2[^ 3) ^ (2) ]-s t2 - tl [0 (2) ^ (1) ] (4) 



where the configurations <f>W specify states in the field diagonal representation at time tj and the exp(—St-t> [</>, </>']) 
denotes the matrix element of the Euclidean time evolution operator during the time interval t — t'. This equation can 
obviously be derived for any theory with nearest neighbor coupling in time. New variables <j>(x) — > 4> a (x), a = 1, . . . , rid 
which allow us to rewrite the action with higher derivative in a form with nearest neighbor coupling in time can be 
introduced in the following manner. Start with a hyper-cubic lattice with lattice spacing a = 1 in each direction and 
construct an anisotropic lattice where the lattice spacing in the time direction is increased to rid by regrouping rid 
time slices of the original lattice. A natural choice is (f) a (x) = <9q0(x), the a-th order finite difference operator in time 
acting on the original field where the finite difference is calculated from the center of the blocked time slice in a time 
reversal covariant manner assuming odd rid. The map <p(x) — > <p a (x) of the Euclidean field variables is an invertible 
linear transformation which preserves the lattice regulated action, Se [4>] — Se \<j> a ] and the generator functional, 

Z E [j] = J D^e-^M+Z 

= II / D[(b a ]e- SE[4 ' a]+ ^ 3 ' t ' (5) 

a J 

as long as the source is placed at the center of the block time slices. The transformation preserves its form in 
Minkowski space-time and provides the mapping whose inverse can be used after the Wick rotation of the blocked 
time slice theory to real time. 

The signature of the norm of states created by the operator (j) a {x) turns out to be <r[<f) a } = (— l) a . In order to preserve 
the orthogonality of field eigenvectors, (4>\4>') = for 0(x) ^ </>'(x) we have to use skew-adjoint field operators, which 
possess imaginary eigenvalues, in the negative norm sector and a\<p\ = ±1 for self- and skew-adjoint variables. It is 
useful to introduce fields with well defined time reversal parity, T(f>(t) — r{(f>](f>(—t) , giving r[d^(f>] = (— l) a r[</>]. This 
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relation suggests the equivalence of the internal Euclidean time reversal parity and the signature of the state created 
by acting on time reversal invariant vacuum by any time reversal invariant combination ip of elementary fields a , 

<#]=t[# (6) 

One has to make sure that unitarity holds within the physical, positive norm subspace, too. This can be achieved 
by the reconstruction theorem of axiomatic quantum field theory, in particular by showing that the main nontrivial 
condition of the theorem, reflection positivity holds in the linear space generated by the action of local operators with 
positive time parity on the vacuum as long as both dynamics and vacuum respect time reversal invariance and the 
boundary conditions 4> a (tf,x) = (— l)°0 a (^,x) are imposed where ti and tf denote the initial and final time. An 
important result of the argument [14] is the direct verification of Eq. (6). This relation indicates, as well, that the 
trajectory of (f> a in the path integral is real or imaginary for a even or odd, respectively. The vacuum may contain 
condensate as long as it is invariant under time reversal. 

This construction gives in the first sight more than expected, it eliminates non-unitarity altogether for theories 
(l)-(2). But the tacit assumptions the argument relies upon are the convergence of the Euclidean path integral 
and the possibility of its analytic extension, Wick rotation, back to real time. The former condition (i) imposes 
> 0. The latter assumption requires that the rotation of the frequency contour in the loop integrals is carried 
out without passing singularities in the integrals. This conditions excludes poles from the quadrant • > 
of the complex energy plane. Since poles come in complex conjugate pairs the remaining complex poles break time 
reversal invariance and generate acausality known from the attempts of removing self acceleration of point charges in 
classical electrodynamics [16]. Thus time reversal invariance restricts the argument to theories where the roots of L 
are real. 

Note that the exclusion of complex poles from the kinetic term restricts the space-time dependence of the pertur- 
bative Green functions to the sum of oscillatory terms e tLOt excluding monotonic terms like e ut . The functional space 
in which the expectation values are constructed is tailored in this manner and the runaway solutions, characteristic of 
unstable theories are excluded. This is in contrast to classical physics where the integration of the equations of motion 
is performed in an unlimited functional space of trajectories. Therefore the classical and the quantum, loop-expansion 
based stability analysis disagree as far as the time-dependent instabilities are concerned. This eliminates the notorious 
instability problem of theories with higher order kinetic term [17]. 



IV. SCALAR ELECTRODYNAMICS 



An important step towards more realistic models is the extension of previous discussion for gauge models. We now 
turn to scalar electrodynamics, defined by the Lagrangian 

C = -\j dxF^F^ + J dx[<j>* L(-D 2 )0 - V(<P*<P)l (7) 

with F^ v = d^By — d u B^ and = — ieB^, L(z) being a polynomial of finite order and is supposed to possess 
separate time and space inversion invariance. In relativistically covariant canonical quantization procedure one adds 
a gauge fixing term, C — > C — £(<9A) 2 /2, and imposes the canonical commutation relations 

[A^t, x), n„(t, y)] = -iff^x - y) (8) 

where IT 1 = dC/ddoA^ and g^ v = (1, — 1, — 1,-1). The wrong sign on the right hand side for [i = v = indicates 
that temporal photon states have negative norm. The Gupta-Bleuler quantization procedure or BRST symmetry can 
be used to prove that usual QED, without higher order derivative is unitary in the physical subspace, spanned by 
states with positive norm. 

With an A field represented by a self-adjoint field operator the field eigenstates are not orthogonal. Orthogonality 
is assured if the operator A is skew-adjoint only [14]. The complication, induced by the use of the traditional self- 
adjoint representation is that non-orthogonality renders the path integral expression for the transition amplitudes 
rather complicated. How to recover then the standard path integral representation for gauge theories in Minkowski 
space-time? The usual path integral over real field configurations A fJ- (x) can easily be found by treating Aq as an 
auxiliary, non-dynamical field either in static temporal or Coulomb gauge. The former will be imposed to establish 
unitarity in the physical subspace because the impact of non-vanishing vacuum expectation value for A on the 
dynamics and the similarity with spontaneous symmetry breaking can better be seen in static temporal gauge. The 
latter gauge will be used in clarifying the physical content of the theory since the dynamical degrees of freedom can 
be traced easier in Coulomb gauge. 
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We start with fields defined without initial or final conditions in time for — oo < t < oo and carry out the gauge 



transformation — > + d^a, cj> — > e lea (f) and 



a(t, x) = - J dt'A (t', x) (9) 

to arrive at temporal gauge Aq = where the functional Schrodinger representation is constructed, using A(x) as 
coordinates. The canonical momentum II = 8qA = — E satisfies the canonical commutation relations [A,(x), Ilfc(y)] = 
iSjkS(x — y). Gauss' law, VE = p, where p is the electric charge density, the equation of motion for Aq is lost in this 
gauge but can be regained as a constraint. In fact, it can easily be shown by the help of the canonical commutation 
relations that 

G[a] = J d 3 i[V«(x)E(x) + a(x)p(x)] (10) 

generates static gauge transformations hence it commutes with the gauge invariant Hamiltonian H, [G(x),H] = 0. 
The average over static gauge transformations, 

V = J L)[a]e i ^' i3;l;[VQ(x)E+a(x)p(x)] (11) 

projects into the subspace satisfying Gauss' law for a given static charge distribution p(x). 

One is usually interested in transition amplitude between gauge invariant states, the latter constructed from a 
gauge-noninvariant representative like a field eigenstate, 

\A,<f>,<l>*) avm = P\A,<f>,<f>*). (12) 

It is enough to insert the projection operator V once only in the matrix element, 

= J D[a](A / ,^ / ,^|e^ d3x [ Va W E+a W' , W]e- itff |A i ,^,^*) (13) 

and one finds the path integral representation 

(A f ,cj> f ,<f>* f \e- itH \A i ,<j> i ,<j>*) sym = J D[A]D[tf>]D[tf>*]e iS ^ A ^ (14) 

where S st [A, (j), (j)*] is the usual action in static temporal gauge, 

d A (x) = 0, (15) 

and tA (x) = a(x) denotes the time- independent integral parameter of the projector. The integration is over config- 
urations A(U,x) = Ai(x), 0(ti,x) = <k(x), 4>*(U,x) = 0*(x), A(t/,x) = A/(x), (t>(t f ,x) = 0/(x), 4>*(t fl x) = 0}(x). 
If the projector is inserted at each time slice of the path integral expression for transition amplitude then the gauge 
invariant action is recovered, 

(A f ,cf> f ,<f>}\e- iAtH P---Pe- iAtH \A i ,cf> i ,<f>*) sym = J D[A]D[4>]D[4T]e iS l A 'W (16) 

AtAo(t, x) playing the role of parameter a(x) in the projector inserted at time t. 

Temporal gauge, used in the Hamiltonian formalism after Eq. (9) is usually not accessible when boundary conditions 
are imposed in time, as done in path integral expressions. Actually the field component A (x) represents a true 
physical variable. We can see this by noting that A (x) cannot be transformed away from the path integral by gauge 
transformation. In fact, setting Aq = instead of integrating over Aq(x) on the right hand side of Eq. (16) removes 
the projector P on the left hand side and the matrix element is changed, (• • -) S2/m — > (■■•). 

A generally applicable gauge choice is static temporal gauge, given by Eq. (15). Whatever gauge we use, the 
Polyakov line 

n( x ) = e- ieI *? dtAo{t ' x) (17) 
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denotes a physical, gauge invariant quantity which prevents us from reaching temporal gauge as soon as some boundary 
conditions arc imposed at the initial and final time. But the integrand of the path integral (f 4) remains unchanged 
under global gauge transformation of the initial or final state by the phase factor 1 = exp2iri, represented by the shift 

Ao(x M A (x) + -^. (18) 

Due to this discrete symmetry the integrand in Eq. (16) does not depend on the space-time independent component, 
A (x) = A and the variable A decouples in the limit tf — U — > oo. Nevertheless, the homogeneous component 
A remains a physical parameter when matrix elements among the vacuum are considered because the vacuum state 
depends on A . In fact, eA acts as a chemical potential and one arrives at grand canonical ensemble where expectation 
values of observables are saturated by the total charge sector of the Fock space which minimizes H — eA J d 3 xp. 



V. SEMI-CLASSICAL VACUUM 



Let us suppose that the model given by Eq. (7) is weakly coupled and saddle point expansion can be used to 
explore its phase structure. The case of global symmetry, e = in the absence of higher order derivative terms 
L(p 2 ) = p 2 , is well known, the model supports homogeneous condensate for appropriately chosen local potential. 
Higher order derivative terms in the action may induce a condensation of particles with non- vanishing momentum, 
an inhomogeneous coherent state and a relativistic "band structure", reminiscent of solid state physics is observed. 
When the interaction with the gauge field is turned on with L(p 2 ) = p 2 then the usual Higgs phase can be found. An 
interesting variant of Higgs mechanism can be generated by the higher order derivatives terms. The point is that the 
partial derivatives are turned into covariant derivatives in the minimal coupling scheme and contain the connection 
term which can induce a nontrivial local potential for the gauge field. The effective interaction, represented by this 
potential may induce a non-vanishing expectation value for the gauge field. We call such a vacuum condensate though 
one should keep in mind that it is actually a coherent state only because our gauge particle, the photon is neutral 
and the Bose-Einstein condensation is not possible. 



A. Condensate 



We follow the strategy of the saddle point approximation and for this end we separate the fields into the sum of 
saddle point and quantum fluctuations by writing </> = 4> + X an d = A^ + A^, the first term in each expression 
representing the saddle point. When a non- vanishing value of the covariant derivative 

-D 2 4>(x) = k 2 4>(x) (19) 

is selected for the semi-classical vacuum by the kinetic energy of the charges then a gauge transformation can always 
exchange contributions of the partial derivative and the connection term. One possibility is when the eigenvalue k 2 in 
this equation is provided by the partial derivative alone, 4>(x) = 4>e~ lpx , A^ = 0. By a suitable gauge transformation 
we may rearrange the semi-classical vacuum into <f>(x) = <j>, eA^ = k^. This is a remarkable simplification offered 
by gauge invariance, the vacuum consisting of the condensate of particles of non- vanishing momentum can be made 
homogeneous. We exploit this possibility and assume the homogeneity of the saddle point and the orthogonality of 
the fluctuations to the saddle point, 

J dxx(x) = J dxA^x) = 0. (20) 

Note that apart from broken global gauge invariance the gauge field condensate leads to the spontaneous breakdown 
of the Lorentz symmetry. When the function L(p 2 ) generates spacelike gauge field condensate, k 2 < then Lorentz 
symmetry is reduced to 0(1,2) and the excitation spectrum loses rotational invariance. We seek vacuum with non- 
relativistic Galilean 0(3) invariance hence we restrict our attention to models with timelike gauge field condensate, 
eA^ = g^ok > 0. Hence there will be four combinations of fields playing the role of Goldstone bosons when <f>, A^ ^ 0, 
corresponding to gauge rotations and Lorentz boosts. The number of massless particle modes is not necessarily the 
same. On one hand, it may be smaller because either non-relativistic fields have half as many particle modes as 
their relativistic counterparts [18] or some field combinations may control not particle-like excitations, with vanishing 
residuum in the propagator at the "mass-shell" . On the other hand it may be more because higher order derivative 
terms may generate several "bands" . Global gauge rotation is applied if necessary to make the scalar condensate, cf>, 
real. 
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B. Fluctuations 



According to section III classical stability analysis is sufficient for the homogeneous components of the fields and the 
stability of the fluctuations around the vacuum will be verified by checking the spectrum of the elementary excitations 
in quantum theory. The energy-momentum tensor of a theory with polynomial, higher order derivative terms can 
easily be obtained, it is the sum of the usual expression for the energy-momentum tensor plus terms containing higher 
order derivatives of the fields. Therefore the energy density of the semi-classical homogeneous vacuum characterized 
is given by the Lagrangian up to a sign, 



by and 



U{e z A 2 , a ) = -^L{e z A z ) + V(0 a ). 



(21) 



We assume at this point that L(p 2 ) is bounded from above and it assumes a maximal value at p 2 = k 2 thus the 
minimization with respect to A 2 , 



n dU(e 2 A 2 ,cj> 2 ) 2 J2N 

° = de 2 A 2 <t>L{eA) 



(22) 



sets e 2 A 2 = k 2 and eA^ — g^fik as mentioned above. The separation of the kinetic and the potential energy 
term in the Lagrangian (7) for the scalar field is not unique, the invariance of the action under the transformation 
L(p 2 ) -> L{p 2 ) + AL, V{4> 2 ) -> V{<p 2 ) - AL<f> 2 can be used to set L(k 2 ) = 0. We assume the form 



L(p 2 ) = -^( P 2 -k 2 ) 2 , 



(23) 



the simplest polynomial satisfying our requirements. The scalar condensate <f) is found by minimizing U(k 2 ,(f> 2 ), i.e. 
solving the equation 



= V'{4> 2 ) - L(k 2 ) 



(24) 



with the auxiliary condition that the first non-vanishing derivative of the potential at the vacuum is positive. 

Once the homogeneous field components are found we turn to the free theory by considering the quadratic part of 
the action. We use the decomposition \ — Xi + *X2, and A = nA^ + At, n = p/|p|, followed by the separation of 
the static components \ a , Al and At by writing Xa — > Xa + Xa, Al — > Al+ Al and At — > Ay + Ay. The quadratic 
action is written as a sum = + with 



S (2) = \f d i x(xu X 2,A L ,A T ) 





K 12 


K 1L 
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\A t ) 
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The momentum space representation of the quadratic form, 

K{p) = J dxe ipi - x - y) K{x,y) 

is calculated in Appendix A with the result 



Ku 


= L+(p)-AV"4> 2 -- 


= i^22 


K 12 


= iL d(P) = ~ K 21 




K 1L 


= -\p\[z(p)L d (p)}- 


= K L1 


K 2L 


= i\p\[z(p)L d (p)] + 


= -K L2 


K LL 


= p 2 [z 2 (p)L d (p)} + 




Ktt 


= - 2 - P 2 , 





(25) 



(26) 



(27) 
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where the notation f ± (p) = f(p) ± f(—p) has been introduced with 

L d (p)=L((p + eA) 2 )-L(k 2 ), (28) 

z{p) = ecf>/ {p 2 + 2ujk) and p = (w, p) for the four dimensional fields. The three dimensional, static sector has quadratic 
forms K(p) = K(p)\ u=0 , obtained from Eqs. (27) and 

K> 4e< ^ fc T f \ K- 

K w = —L d {p)\ lJJ=0 = K 0l 

= ( p 2)2 Mp)|^o + P 2 (29) 



C. Unitarity 

We turn now to the question of unitarity of the time evolution within the positive norm subspace of the Fock-space. 
There are two circumstances requiring to go beyond the argument based on the reconstruction theorem for Euclidean 
theories [19]. One is that the manifest 0(4)/Lorentz invariance of the Euclidean/Minkowski Green functions, one 
of the numerous conditions of the theorem is lost in our case. Another point is that states belonging to excitations 
generated by the time component of the gauge field have negative norm in Minkowski space-time and are thus non- 
physical. Rather than attempting to generalize the reconstruction theorem we choose a simpler argument, valid in 
any finite order of the perturbation expansion. 

The partial fraction decomposition of the propagator is now made in terms of lu 2 rather than p 2 and the realness 
of the one-particle energies guarantees the unitarity of the perturbative model within the Fock-space with indefinite 
norm. Perturbation expansion, based on the vacuum with homogeneous fields (f> and leads to a stable and unitary 
theory if all solutions of the equation det K(j>) = of the quadratic form K of Eq. (25), 



det K(p) 



fc 4 



2 p 2 ) 2 ^ 10 



{• 



4w 8 (p 2 + 2fc 2 ) + lu 6 [l6fc 4 + 16fc 2 p 2 + 6 (p 2 ) 2 + W"^ 2 k 2 



W"0 2 (e 2 p 2 2 + 2fc 2 p 2 - 4fc 4 ) + 4 (p 2 ) 3 + 8k 2 (p 2 ) 



6W" 2 4 + W"4> 2 k 2 (p 2 ) z + 2e 2 4> 2 (p 2 ) z - 4e 2 2 fc 2 p 5 
-2e 2 k 2 p 2 V"0 3 - 4e 2 (p 2 ) V> 3 } , 



(30) 



obtained for the kinetic term (23) have real frequency components, lo 2 > 0. It is easy to see that this expression has 
negative or complex lu 2 as roots, there are instable modes in the scalar particle, longitudinal gauge field sector. These 
instabilities can be excluded by imposing the condition 







(31) 



on the local potential which is not a natural relation, it requires a fine tuning to cancel the scalar particle scattering 
amplitude at vanishing momentum. 

According to Goldstone theorem the minimization of the vacuum energy with respect to the strength of condensate 
cancels the gap for certain modes. Goldstone mode arising from the breakdown of global gauge invariance is made by 
Eq. (24). As far as the three soft field combinations are concerned, which correspond to the breakdown of Lorentz 
symmetry, let us introduce a mass term for the gauge field by the extension C — >• C + m 2 B 2 /2 of the Lagrangian (7) 
as in Proca theory which leads to the modified potential U(e 2 A 2 ,(j) 2 ) — > U(e 2 A 2 ,<f> 2 ) — m 2 A 2 /2 in Eq. (21). The 
minimization with respect to the gauge field condensate, Eq. (22), generates three soft field combinations. Two of 
them are the non- vanishing helicity components of transverse gauge field even for m 2 ^ and a third is a combination 
of dpA* 1 , A^A^, xi and \i- To simplify matters we return in our discussion to scalar electrodynamics, m 2 — where 
the determinant of Eq. (30), whose vanishing identifies the normal mode dispersion relation, reads 



A C tK{p) = Yi {u 2 



,2\2, ,2/ 



w 2 (w 2 - 2kcu - p 2 ) 2 (u; 2 + 2ku - p 2 ) 2 . 



(32) 



The energy spectrum is real, transverse gauge fields make up two Goldstone modes with oj = ±|p|. The scalar field, 
together with the longitudinal components of the gauge field produce the dispersion relations 



u\k + <r 2 \/k 2 + p 2 



(33) 
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where cti,CT2 = ±1. The choice <j\(J2 = — 1 in Eq. (33) belongs to two other Goldstone modes. The determinant (30) 
correspond to non-static fluctuations, therefore the factor uj — in Eq. (32) is never vanishing. 

Once the unitarity has been established in the whole Fock-space let us turn to the physical subspace. The argument 
in Ref. [14] was presented for Yang-Mills-Higgs model, given by the Lagrangian (7) though some additional care is 
required in this case to draw conclusions for Minkowski space-time theories. The Wick rotation is more involved 
for gauge than for scalar fields because the norm of state created by Ao changes sign during Wick rotation between 
Euclidean and Minkowski space-time. This leads to the following two problems. One has already been mentioned in 
Section IV, the usual path integration formulas require the orthogonality of the field eigenstates and we should use 
skew-adjoint representation for Aq in Minkowski space-time. This amounts to integration over imaginary Aq field 
which is in an obvious conflict with the usual interpretation of A as the temporal component of a Hcrmitian quantum 
field. The solution of this apparent contradiction is well known, the treatment of A as a non-dynamical, auxiliary 
variable. This is what happens in static temporal gauge where Aq is the (real) integral variable of the projection 
operator to restrict the dynamics into the subspace with Gauss' law. Once the real, static Ao configurations are 
accepted in the path integral of Eq. (14) then we may return to the gauge-free case, Eq. (16) in the calculation 
of gauge invariant quantities. In other words, in the usual path integral formalism for real time, available for gauge 
theories with higher order derivatives, as well, the temporal component of the gauge field is better to interpret as 
an auxilary variable to handle Gauss' law rather than a quantum field handling physical excitations. The situation 
is reminiscent of conventional QED where elementary excitations, stability, renormalizability, etc. are trivial in 
relativistic gauges but one has to go into another, physical gauge, usually chosen to be the Coulomb gauge to recover 
unitarity in the physical subspace in an obvious manner. 

Another problem, caused by an exceptional feature of A (x) during Wick rotation is that Eq. (6) used to identify 
the signature of the norm is not valid anymore for this component of the gauge field in Minkowski space-time. 
A generalization valid for gauge field is 



where ip is any local combination of the elementary bosonic fields d§<j>* and dy A and space inversion acts as 
Pip(t,x) — ip(t, — x) with tt[4>] — — 7r[A] = 1. PT invariance yields the conservation of a and assures unitarity within 
the positive norm, physical subspace [21]. 



It has been established so far that our model has unitary time evolution within the positive norm subspace and is 
therefore a physically interpretable. The next question is its physical content which will be assessed by comparing it 
with standard electrodynamics. The usual Higgs-mechanism renders photons massive. The Goldstone modes arising 
from the spontaneous breakdown of the Lorentz invariance make three combinations of the fields soft. Two of them 
are the transverse, non-vanishing helicity components of the gauge field and they keep the radiation field massless, 
just as in standard electrodynamics. Two further soft field combinations are made up from the longitudinal gauge 
and the scalar field components. 

The double pole of (23) may render the normal modes of the scalar field non particle-like because scattering 
amplitude wave packets, constructed by this kind of excitations may be vanishing according to the reduction formulas. 
Thus we take the point of view that the scalar field corresponds to so far non-observed excitations and seek the 
dynamics of the gauge field only. To simplify matters further, we ignore radiative corrections due to the charged 
scalar field and restrict ourselves to the O (A 2 ) part of the action where the normal modes are quasi-photons. Wc 
consider below two aspects of the model, the number of propagating, dynamical degrees of freedom and their dispersion 
relation. It is worthwhile separating two different kinds of dynamics for the gauge field, first arising through the field 
strength tensor in the Maxwell-action, the first term in the Lagrangian (7) and second, coming directly from the 
connection term of the covariant derivative in the minimal coupling. The former, field strength tensor dynamics 
represents conventional electrodynamics and the latter, connection term dynamics is the source of genuine quantum 
and topological effects. 

Let us first have a look into the Proca-theory, the simplest model with massive vector field and use the standard 
three-dimensional notation A^ = (tp,A), — (p,j), E = — Vip — 8qA, H = V x A. We separate the transverse 
and longitudinal components, A = At + V$, j = jt + Vk where current conservation implies dop + Ak = 0. The 
Lagrangian 



a[lp] — T[l/>]7r[?/>] 



(34) 



VI. QUASI-PHOTONS 



£ = ±E 2 -iB 2 - w +^ 2 -A 2 )+jA 



(35) 
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can be written as L — Lt + Llq where the first and the second term contain the transverse and longitudinal and 
temporal components, 



I TYI, 

£lo = -(V^ + a V$) 2 + — [<^ 2 -(V$) 2 ]- W + V$Vk. (36) 

As is well known, the temporal component ip is not a dynamical degree of freedom and can be eliminated by solving 
its algebraic equation of motion in time, 

V=^—r(p + d a A3>), (37) 
m z — A 

without generating non-local effects in time and the resulting effective Lagrangian for $ is 



1 1 to 2 A(D + m 2 l ^ 

2 m z — A 2 ra z — A 



A 



-d p — An 



(38) 



For massless photon, m 2 = 0, the equation of motion for $ is the current conservation and longitudinal photons drop 
out from the field strength dynamics. But the mass term, arising from the connection term dynamics may bring the 
longitudinal component back as a genuine dynamical variable. Gauge transformations may make the separation of 
auxiliary and truly dynamical variables difficult. For instance, there are gauges, such as the static temporal gauge, 
where the longitudinal component appears to be dynamical but it drops out from gauge invariant observables. When 
higher order derivative terms appear in the connection term dynamics then either the temporal or the transverse 
component of the gauge field may acquire non-trivial dynamics. The formal gauge invariance always makes the theory 
redundant therefore one expects three dynamical, propagating degrees of freedom for the theory (7) from the photon 
field, just as in the usual Higgs-mechanism. But their dispersion relations differ from those of the Higgs-mechanism, 
betraying the different underlying symmetry breaking patterns. 

Let us look into the dispersion relation of the model (7) in Coulomb gauge which offers a particularly clear view in 
our model with spontaneously broken Lorentz symmetry. The Lagrangian L = L T + L 0m is written as the sum of the 
transverse part, given by the first equation in Eqs. (36), and the rest whose quadratic part is 



4m = n(Xi,X2, V)K C | X2 ) (39) 



where X — Xi + *X2 and 



L+(p) iL-{p) [(p° + 2k)z(p)L d (p)]+ 

K c =\ -iL- d {p) L+(p) -i[(p° + 2k)z(p)L d (p)]- ). (40) 

K [(p° + 2k)z(p)L d (p)}+ l [(p" + 2k)z(p)L d (p)}- [{2k+p*Yz 2 {p)L d (p)]+ +v\ 

The dispersion relation is defined by the roots of the determinant of the quadratic form, 

det K c {p) = jJipV 2 - 2fc^ - P 2 ) V 2 + 2kco - p 2 ) 2 . (41) 

Comparing this expression with Eq. (32), the determinant of the small fluctuations in static temporal gauge apart 
from the obvious absence of two massless modes, corresponding to non- vanishing helicity transverse modes of the gauge 
field one notices the appearance of a new root, p 2 , suggesting the emergence of a conventional Coulomb propagator. 
One can obtain a more detailed view of the normal modes by the inspection of the propagators. The inverse of Kc 
is a full matrix with rather involved matrix elements. Matrix elements of Kq 1 between the matter field contain the 
factor (cj 2 — 2fcw — p 2 ) 2 (w 2 + 2kiv — p 2 ) 2 , indicating the non-particle like behavior. The matrix elements between the 
matter field and tp have the factors (co 2 — 2ku) — p 2 )(w 2 + 2ku> — p 2 ) and p 2 in the denominators. Finally, the simplest 
inverse matrix element is the diagonal one for ip, 



(Kc'U = 4j, (42) 



P 2 



confirming that the factor p 2 in Eq. (41) corresponds to the unchanged Coulomb law. 
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The expectation of three dynamical, propagating components for the gauge field, mentioned after Eq. (38) above 
turned out to be wrong and the nontrivial dynamics for the longitudinal component, expected by analogy with the 
Proca case, Eq. (38) was too naive. The higher order derivative terms render the nontrivial dispersion relation for the 
longitudinal component a gauge artifact and the usual dispersion relation is recovered for the electromagnetic field. 

The surprising simplicity of (42) is the result of nontrivial cancellations. This can easiest be seen by calculating 
the Av propagator directly. For this end we eliminate the charged field by its equation of motion what is simplest to 
carry out in the complex x basis, where 

1 (K- K \ f X \ 

C^ = -{x\x^)\ X * (43) 

2 \K K oa J \tpj 

with K_ = -2(D-2ikd a ) 2 /k 2 , K = 2ej){2k + id ){U - 2ikd ) / k 2 and L 00 = 2e 2 4> 2 (d% - 4k 2 )/k 2 - A. The equations 
of motion for x* and X, = K_\ + K a A a , and = x^K- + ^0-^0, used to eliminate the scalar field yield 

£^ = \<pD^<p (44) 

where 

= K m - ^(KoKZ'Ko + Kl r KZ ltr K) (45) 

gives Dq Q = p 2 after some cancellations. Therefore the deviation from usual electrodynamics and the impact of the 
higher order derivative terms are seen by the charged scalar field in our approximation. 



VII. CONCLUSION 



A novel spontaneous symmetry breaking is discussed in the framework of scalar QED which involves higher order 
covariant derivatives. One finds non- vanishing expectation value for the gauge field and unitary, physically acceptable 
interactions in properly fine tuned models. 

The unitarity is proven in the physical, positive norm subspace in two steps. First it is assured in the whole Fock 
space by fine tuning the self interactions for the charged scalar field. Second, it is shown that PT invariance makes 
the physical subspace closed under time evolution. 

The particle content of our model is radically different than the one found in the conventional Higgs-mechanism. 
Goldstone theorem renders the radiation field massless. Furthermore, a particular model is proposed where all 
components of the gauge field are massless and Maxwell equations are recovered in the linearized equations of motion. 

We sought in this work a vacuum which supports Galilean invariance, therefore the temporal component of the 
gauge field was allowed to develop vacuum expectation value. It acts as some dynamically generated chemical potential 
for the charged scalar particle. The scalar particle condensate remains electrically neutral due to the equal number 
of particles and anti-particles it contains as a result of the higher order derivative terms in their dispersion relation. 

The status of Lorentz symmetry, broken by the vacuum expectation values to the Galilean group, is rather peculiar 
in the Abelian model. Despite the breakdown of Lorentz invariance Goldstone modes display relativistic dispersion 
relations. Furthermore, three components of the gauge field become Goldstone modes corresponding to the spontenous 
breakdown of relativistic symmetries and hence remain massless even if one starts with massive Proca action for 
photons. The quadratic part of the Lagrangian in the fluctuations of the gauge field is identical to that of QED, 
leaving the Lorentz non-invariant part of the photon dynamics to be generated by radiative corrections. The deviation 
of this model from standard electrodynamics is due to radiative corrections only. 

There are numerous extensions one may consider. Similar models with non-Abelian gauge symmetry should lead 
to some massive gauge field components because Goldstone theorem cannot protect all components of the gauge field 
anymore against mass generation. Using a basis in internal space where the massless gauge bosons are diagonal the 
other, non-commuting components of the gauge field are charged and allow us to construct models with an unbroken 
U{\) subgroup, as in the Standard Model. It remains to be seen if natural models, requiring no fine-tuning can be 
constructed by the eventual inclusion of charged fermions. Another issue, the scale-dependence of the breakdown of 
Lorentz invariance is interesting, too. Being a spontaneous symmetry breaking, it should be strong at low energy. 
But some interesting results about non-Lorentz invariant quadratic terms in gauge theories [20] suggests that certain 
Lorentz symmetry breaking parameters of the dynamics tend to be suppressed in the low energy limit. A systematic 
renormalization group study of the model would be needed to reveal the true scale dependence of this symmetry 
breaking. Finally, extension for gravity opens new questions since the spontaneous breakdown of Lorentz symmetry 
may generate massive gravitons by a gravitational Higgs-effect. 
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Appendix A: Quadratic action in momentum space 



To find the momentum dependence of the quadratic form K(p) we evaluate the quadratic action (25) for the test 
functions 

X (x) = X 'e~ ipx , 

A^x) = A'^e-** (Al) 

before gauge fixing for the sake of simplicity. The quadratic form of the O (x*X) P ar t can easily be written as 

K x , x = 2L d (p) - W"4> 2 (A2) 

by means of Eq. (24) with Ld(p) introduced in Eq. (28). 

To find the other terms it is advantageous to represent the higher derivative kinetic term of the scalar field as a 
polynomial, 



„2n 



(A3) 



L (p 2 ) = ° n P 

The block that mixes the scalar and the gauge field originates from the O (B) piece in 

rid 

L{-D 2 )4> = J2 c n[-(d - leA - ieA) 2 ] n 4> 
and we find for the O (A\) contributions 

- / dxdy X (x)K xA (x, y)A{y) = ie dx X (x) £ c n • ■ • (2A(x)B + 8A{x)) ■ ■ ■ (-0)0 (A5) 



(A4) 



n=0 e=i 



where = <9 - ieA^, □ = d^, and the l-th factor of the term O ((-D 2 ) n ) is replaced by the O (A) part of -D 2 
in the right hand side. The choice (Al) leads to 



written as 



K xA {p)A' = 2ie °n + eA ) 2 ' ' ' (- 2iA o k - ipA')k 2 

n=0 1=1 



K xA (p)A> = 2e(2A> k+pA')k- 2 £ c n k 2n £ 1 + P + f 

n=0 t=n ^ ' ' 



The geometric series can be summed, 



p 2 + 2p°k 



n=0 



and we have 



K xAfi (p)=2e0 



L d {p) 



2g»°k + pi" 



p 2 + 2p°k 



(A6) 



(A7) 



(A8) 



(A9) 



The O (A 2 ) quadratic form for real field requires more care. Since it acts on real field it must be symmetrical. We 
shall consider a complex plane wave component of the gauge field in the actual calculation of this term and carry out 
the symmetrization at the end only. This term is the sum of two contributions. One of them is the standard Maxwell 
piece 



v m 2 



(AlO) 
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where = g^ v — p^p v jp 2 is the projection into the transverse polarization subspace. The other part is the O {A 2 ) 
contribution in 



4>L{-D 2 )4> = Cn<f>[~(d - ieA - ieA) 2 \ 



(All) 



which will be written as the sum K AA (p) + K AA (p). The first term stands for the O (A 2 ) contributions of the — D 2 
factor, 



A'Kf A A> = 2e 2 ^ £ #-□') • • • A 2 {x) ■ ■ ■ (-□> 



n=0 e=i 



which is vanishing, 



KW(p)=2g>»'4?e 2 L'(k?) = 0. 



what is written as 



rid 7i—l n 

A'K A 3 A (p)A' = 2e 2 $ 2 (2A k+ P A) 2 k- i ^Tc n k 2n J2 E ( 1 + 

1=1 e'=e+i 



n=0 



p 2 + 2p a k 
k 2 



e'-e-i 



The summation of this geometric series gives 

,(2A k + P A) 2 



n=0 



n-1 



p 2 + 2p° k" n 
k 2 



n + 1 



The resulting geometrical series in the square bracket can again be summed with the result 



1 / fi '(3)^>,/_o„2^(2A)fc+pA) 2 



A'K%Sp)A' = 2e 2 



yielding finally 



and 



(p 2 + 2p°k) 2 



A'K A 3 A (p)A' = 2e 2 4> 2 L d {p) 



^ ' c n k 



n=0 



p 2 + 2p°k 

k 2 



{2A k+pAf 
(p 2 + 2p°k) 2 



l-n 



p 2 + 2p°k 
"fc 2 



K A ^(p) = -T^ + e 2 4> 2 L d (p) 
after symmctrization. 



• - V e (p L d {-p) — 



(p 2 + 2p°k) 2 



(p 2 - 2p°k) 2 



(A12) 



(A13) 



The other contributions is for the product of two O (A) terms, 

n— 1 n 

A'Kf A {p)A' = -2e 2 J2c n J2 '-O • • • (2Ad' + dA) ■ ■ ■ (2Ad' + dA) ■ ■ ■ (—□')(/> 

n=0 £=1 £'=£+1 
rid n—1 n 

= 2e 2 ^c„^ 4>k 2 ■■■(2A k+pA)---(k 2 +p 2 + 2p°k)---(2A k+pA)---k 2 4> (A14) 

n=0 e=i l>=l+l 



(A15) 



(A16) 



(A17) 



(A18) 



(A19) 
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